Pressure-dependence of electron-phonon coupling and the superconducting phase in 

hep Fe - a hnear response study 
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A recent experiment by Shimizu et al. (Ref.^) has provided evidence of a superconducting phase 
in hep Fe under pressure. To study the pressure-dependence of this superconducting phase we have 
calculated the phonon frequencies and the electron-phonon coupling in hep Fe as a function of the 
lattice parameter, using the linear response (LR) scheme and the full potential linear mufhn-tin 
orbital (FP-LMTO) method. Calculated phonon spectra and the Eliashberg functions a^F indicate 
that conventional s-wave electron-phonon coupling can definitely account for the appearance of the 
superconducting phase in hep Fe. However, the observed change in the transition temperature with 
increasing pressure is far too rapid compared with the calculated results. For comparison with 
the linear response results, we have computed the electron-phonon coupling also by using the rigid 
mufBn-tin (RMT) approximation. From both the LR and the RMT results it appears that electron- 
phonon interaction alone cannot explain the small range of volume over which superconductivity is 
observed. It is shown that ferromagnetic/antiferromagnetic spin fluctuations as well as scattering 
from magnetic impurities (spin-ordered clusters) can account for the observed values of the transition 
temperatures but cannot substantially improve the agreeemnt between the calculated and observed 
presure/ volume range of the superconducting phase. A simplified treatment of p-wave pairing leads 
to extremely small (< lO"'^ K) transition temperatures. Thus our calculations seem to rule out 
both s- and p- wave superconductivity in hep Fe. 

PACS numbers: 74.70.Ad, 71.20.Bc, 74.20.Mn, 74.90.-|-n 
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I. INTRODUCTION 

Recently Shimizu et aM (see also Ref.||) have reported 
resistivity and magnetization measurements on Fe sam- 
ples under pressure, and identified a superconducting 
phase characterized by the Meissner effect and the van- 
ishing of the resistivity above a pressure of 15 GPa. At 
this pressure the stable crystal sructure of Fe is known 
to be hep. Both the hep phase and superconductivity 
in Fe under pressure are results that can be expected on 
theoretical grounds. Stability of bcc, fee and hep crys- 
tal structures as a function of canonical |d-band filling 
was discussed some time back by PettiforEl and Ander- 
sen et al.nB These authors showed that without ferro- 
magnetism the ground state of Fe would be hep, just as 
for its nonmagnetic and isoclectronic Ad and 5d counter- 
parts, Ru and Os. For Fe the bcc structure is stabilized 
only via ferromagnetism. In the ferromagnetic bcc state 
both the atomic Vjolume and compressibility of Fe are 
anomalously large. □ Application of a moderate pressure 
results in bcc to hcxi martensitic transformation and loss 
of fcrromagnetism.Q Both Ru and Os arc superconduct- 
ing at low temperatures. Thus superconductivity in hep 
Fe is hardly surprising. 

What differentiates hep Fe from Ru and Os is the pres- 
ence of spin fiuctuations. Both ferromagnetic and anti- 
ferromagnetic spin fluctuations are known to suppress 
superconductivity mediated via s-wave electron-phonon 
coupling. A notable example, where ferromagnetic spin 
fluctuations (paramagnons) are believed to suppress su- 
perconductivity completely, is fee Pd. A large density of 
states (DOS) at the Fermi level, A''(0), in fee Pd causes a 
large Stoner-enhanced paramagnetic susceptibility, lead- 



ing to strong ferromagnetic spin fluctuations. Disorder- 
induced superconductivity in fee Pd, due mainly to the 
reduction in iV(0) and therefore in spin fluctuations, 
has been claiiaed experimentally as well as discussed 
theoretically .Erffl A similar effect could conceivably be 
achieved in fee Pd under pressure, but is yet to be ob- 
served. The case for hep Fe is somewhat different, since it 
is believftd,to be close to antiferromagnetidij or cpHiplex 
magneticrJ instability. It was noted by Wohlfarthli3 that 
at the lowest pressures (~10 GPa) at which hep Fe is sta- 
ble, it should be close to an antiferromagnetie instability. 
He also suggested that the antiferromagnetie spin fluctu- 
ations might not be strong enough to suppress supercon- 
ductivity in hep Fe, particularly at elevated pressures, 
where reduction in iV(0) would cause spin fluctuations 
to eventually disappear. Antiferromagnetie spin fluctua- 
tions suppress s-wave superconductivity, while contribut- 
ing to p/d-wave superconductivity. At present experi- 
mental evidence regarding the type of superconductivity 
(s-wave or otherwise) in hep Fe is lacking. 

One can estimate the Tc in hep Fe by using simple seal- 
ing arguments and the observed superconducting transi- 
tion temperature Tc of Ru (0.5 K) or Os (0.7 K) at normal 
pressure. Let us ignore spin fluctuations and consider the 
McMillan expression: 
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where O/j is the Debye temperature, fi* is the Coulomb 
pseudopotential, and Xph is the electron-phonon coupling 
constant, given by Xph — N [0) {P) / M {lo"^) . Consider- 
ing /i* = 0.1, we get Xph = 0.32 for Ru (O^ = 600 
K, see Refj6|). To estimate Xph for hep Fe, we assume 
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that the mean square electron-phonon(ion) matrix ele- 
ment (/^) and the effective spring constant M{iu'^) are 
nearly the same as in hep Ru. The average phonon fre- 
quency, and thus Qdi should then scale as the inverse 
square root of the ratio of the atomic masses, and Xph 
should scale according to the ratio of N{0). The quan- 
tity N(0) can be easily calculated for elemental solids. 
However, we only need to estimate the ratio of this quan- 
tity between Ru and Fe. We can start by assuming that 
A'^(O) is proportional to the inverse d-band width, which 
can be estimated from the potential parameters of the 
LMTO-ASA (atomic sphere approximation) method.El 
Here the band width is proportinal to the potential pa- 
rameter A. Both A jaiid its volume derivative have al- 
ready been tabulatedlj for a large number of elemental 
solids. From the c?-orbital values of the parameter A, 
N{0)Fe/N{0)R., ~ ArJAfc = 539/280 = 1.925. This 
would give Xph =0.62 for Fe, resulting in a transition tem- 
perature of ~ 17 K for hep Fe. If we use published values 
of N{0), for hep FeEJ (corresponding to a pressure P ^ 
10 GPa) and for RrtJ (corresponding to normal prssure), 
then N{0)Fe/N{0)Ru = 20.8/11.8 = 1.76, and Xph=0.56 
for hep Fe. This ratio yields an improved value Tc{Fe) 
= 12 K. Using the measured value of Qd iur-iicp Fe (~ 
500 K at - 10 GPa, see Ref.|l|) gives a valued of 7.6 K. 

For a quick estimate of the pressure-dependence of Tc 
we use a simplified version of Eq.|l|: 



According to the study by Shimizu 
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and resort to the tabulated values of the logarithmic 
derivative of the potential parameter A with respect to 
atomic sphere radius, s (Ref.p^). Neglecting the vol- 
ume(pressure) dependence of the quantities (/^) and fj,* 
in Eqj^, we obtain, for the logarithmic derivative of Tc 
with respect to the system volume V: 
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where jc is the Griineisen parameter. We have used the 
approximations: 
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With the assumption A^(0) ~ A ^, where A is 
the d-orbital band width parameter in LMTO-ASA, 
d\nN{0)/d\nV = -(l/3)dln A/dln s. For the d-orhitals 
of Fe din A/dlns = -4.6. From the reported valuet3 of 
7G = 1.5 in hep Fe, we obtain a value = 6.6 for 

hep Ee, The zero pressure bulk modulus in hep Fe is 165 
GPa.113 The initial (low pressure) logarithmic derivative 
of Tc in hep Fe should thus be close to -6.6/165 (GPa)~^ 
= - 4% /GPa. 

Exercises such as the one outlined above are useful in 
obtaining order of magnitude estimates and in under- 
standing the trend from one element to the next. How- 
ever, quantitative agreement with experimental results 



migh^be missing. 

superconductivity in hep Fe appears at around 



et alE 

15 GPa, slightly above the pressure at which the bcc- 
hcp transition takes place. The transition temperature 
grows slowly from below 1 K to about 2 K at ~22 GPa 
and then decreases steaidily, with superconductivity van- 
ishing beyond 30 GPaJla The rate of decrease of Tc is 
too rapid compared with the estimate derived above. In 
order to reproduce the initial increase of Tc with pres- 
sure, as observed in the experiment, it would be neces- 
sary to include the spin fluctuation effects and possible 
volume dependence of the matrix element (/^). With 
inceasing pressure, spin fluctuations are expected to di- 
minish, causing Tc to rise. The electron phonon matrix 
element may also increase with pressure, as the nearest 
neighbor distances become shorter. It would thus be of 
interest to examine to what extent the observed results 
can be explained via a rigorous ab initio calculation. To 
this end we have used the full-potential linear muffin-tin 
orbitals linear reaopase (FP-LMTO-LR) scheme devel- 
oped by SavrasovEjo to calculate the phonon frequencies 
and the electron-phonon coupling in hep Fcpas a func- 
tion of pressure. The Eliashberg equationsj2l in their 
isotropic Fermi surface averaged form, are used to study 
the pressure-dependence of the transition temperature 
Tc, and the superconducting gap A. Effects of both fer- 
romagnetic and antiferromagnetic spin fluctuations and 
the effects of scattering from magnetic impurities are ex- 
plored to accommodate the experimental data as best as 
possible. We also present a simplified treatment of p- 
wave pairing in hep Fe. So far two other theoretical cal- 
culations, related to superconductivit }r . in , hep Fe and its 
pressure dependence, have axm£!ared.EjE3 Our work dif- 
fers from these publicationsOtJ in asmuch as it presents 
a more rigorous first-principles calculation of the phonons 
and the electron-phonon coupling as a function of the lat- 
tice parameter in hep Fe. 



II. ELECTRONIC STRUCTURE 

There is considerable experimental evidence that at 
room temperature the martensitic transition from the 
bee to the hep .phase in iron takes place at a pressure of 
10-15 GPaSiyfel Recently Ekman et aM have studied 
this phase transition using the full-potential linear aug- 
mented plane wave (FP-LAPW) total energy method. 
Their study indicates a first order ferromagnetic bee to 
nonmagnetic hep transition. These authors carried out 
spin-polarized density functional calculations using the 
generalized-gradient— approximation (GGA) qt, Perdew 
and Wang (GGAl).EZl Steinle-Neumann et al.jtB using a 
later version oLthe GGA by Perdew, Burke, and Ernzer- 
hofer (GGA2),Ea find an antiferromagnetic ground state 
for hep Fe and show that this version of the GGA better 
reproduces the observed elastic properties of hep Fe un- 
der pressure. The possibility of noncoUinear magnetism 
in hep Fe below 50 GPa has also been suggested^ These 
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results are at variance with the Mossbauer study of hep 
Fe under pressure ,-4^ ch has failed to reveal any local 
magnetic moment .ESEa The possibility remains that hep 
Fe stays close to antiferromagnetic or complex magnetic 
instability. In this work we assume a nonmagnetic phase 
for hep Fe under pressurCp-and present results that were 
obtained by using GGAlB Our calculations, using var- 
ious forms of GGA, show that the nature of the ground 
state at various lattice parameters depends very much 
on the exchange correlation potential and the c/a ratio. 
In addition, the energy differences between nonmagnetic, 
ferromagnetic, apd antiferromagnetic (AFM I and AFM 
II configurationsliJ ) states are often small (almost within 
the errors of the method) , as well as dependent on tech- 
nical details, such as the number of k-points in the irre- 
ducible Brillouin zone and method of BZ integration (see 
further discussion in section V). 

According to Ekman et aZ.c3 the bee to hep transition 
leads to a phase with a c/a ratio of 1.57. Our LMTO- 
ASA calculations yield a smaller (by ^ 0.8 mRy) hep 
ground state energy for c/a = 1.57 than for the ideal 
close packing value. In the FP-LMTO calculations the 
difference in the ground state energies for the two c/a val- 
ues is smaller than 0.2 mRy. Previous theoretical stud- 
ies for iron indicate a ver}; spall dependence of the to- 
tal energy on the c/a raticE3£d and in addition, the c/a 
ratio is likely to change with pressure. We have thus 
adopted the simplest option, as in Ref.|l^, and carried 
out all our calculations with the ideal close-packing case: 
c/a = ^y8/3. The electronic stcttcture was calculated 
using Savrasov's FP-LMTO coda23 with a triple-K spd 
LMTO basis for the valence bands. 3s— and 3p-semicore 
states were treated as valence states in separate energy 
windows. The charge densities and potentials were rep- 
resented by spherical harmonics with I < 6 inside the 
nonoverlapping MT spheres and by plane waves with en- 
ergies < 141 Ry in the interstitial region. Brillouin zone 
(BZ) integratiopa were performed with the full-cell tetra- 
hedron methodcj using 793 k-points in the irreducible 
zone. Band structures of hep Fe obtained via the FP- 
LMTO method for all lattice parameters considered are 
in good agreement with the LMTO- AS A bands. FP- 
LMTO bands for the ideal c/a ratio and lattice parame- 
ter of 4.6 a.u. are shown in Fig. ^. 

Table | shows the lattice parameters used in our calcu- 
lations together with the atomic volumes and some calcu- 
lated properties. Pressure {—dE/dV) and bulk modulus 
were calculated by fitting the energy vs. lattice param- 
eter curve, jta the generalized Birch-Murnaghan equation 
of stateEj't3 The equilibrium atomic volume and bulk 
modulus, 69.4 a.u. and 290 GPa, compare well with the 
values obtampd by the FP-LAPW calculations of Ek- 
man et alM (68.94 a.u. and 263 GPa). Calculated 
pressure and bulk modulus values become progressively 
less reliable away from the equilibrium volume. In or- 
der to calculate the Stoner parameter / we introduced a 
small splitting in the self-consistent paramagnetic bands 
by adding small up and downward shifts to the band- 
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FIG. 1: FP-LMTO energy bands in hep Fe for the ideal c/a 
value (-^8/3), and a= 4.6 a.u. The horizontal line shows the 
position of the Fermi level, chosen as the zero of energy. 



TABLE L FP-LMTO results for nonmagnetic hep Fe for the 
ideal c/a ratio: a,— lattice parameter (a.u.), Vo = volume 
per atom (a.u.), P=pressure (GPa), S=BuIk Modulus (GPa), 
A^(0)= DOS at the Fermi level (states/(Ry atom spin)), I— 
Stoner parameter (Ry/atom). 



a 


4.0 


4.2 


4.4 


4.5 


4.6 


4.7 


Vo 


45.25 


52.39 


60.23 


64.44 


68.83 


73.41 


P 


350 


162 


56 


26 


2.3 


-14 


B 


1695 


970 


550 


410 


300 


221 


N{Q) 


4.77 


5.79 


7.05 


7.80 


8.59 


9.46 


I 


0.074 


0.074 


0.073 


0.073 


0.073 


0.075 


1/(1 -IN) 


1.54 


1.75 


2.06 


2.32 


2.68 


3.44 



center parameter C in the LMTO-ASA method. After 
making the atom self-consistent the Stoner parameter / 
was calculated from the induced magnetic moment per 
atom /X, assuming proportionality between band-splitting 
and the Stoner parameter: 



/ = Y,Iidi; Ii = {Cj - a^)/fi; Si = Ni{Q)/N{0) , (5) 



where the arrows indicate spin-up and down states and 
iV/(0) and A^(0) are the Z-partial and total DOSs at the 
Fermi level, respectively. This method yields almost the 
same (pressim-independent) value as that obtained by 
Mazin et al.£3 Both our method and that used in Ref.|lJ 
can be called fixed spin-moment method, except that 
Mazin et al. derived / from the second derivative of the 
total energy with respect to the spin- moment. 
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III. LATTICE VIBRATIONS AND 
ELECTRON-PHONON COUPLING 



We used the linear response code of Savrasov '^^1^^ with 
a triple-K LMTO basis set. The dynamical matrix was 
generated for 28 phonon wave vectors in the irreducible 
BZ, corresponding to a mesh of (6,6,6) reciprocal lat- 
tice divisions. The BZ integration for the dynamical 
matrix was done for a mesh of (12,12,12) reciprocal lat- 
tice divisions, and that for the electron-phonon (Hop- 
field) matrix was done for a (24,24,24) mesh. The cal- 
culated phonon spectra for two lattice parameters, 4.4 
a.u. and 4.0 a.u., are shown in Fig.|[ Our results are 
in reasonable agreement adth a recent density functional 
calculation of Alfe et alE3 These authors usje-|the same 
GGA as is used in our c fOculation (GGAl)E!l and the 
small displacement methodEZl to obtain the force constant 
matrix. In Fig. 3 of their paper the phonon spectra for 
two volumes 8.67 A'^ and 6.97 A"^ are shown. The cor- 
responding lattice parameters, 4.36 a.u. and 4.05 a.u., 
are close to the values for which our calculated phonon 
dispersions curves are shown in Fig. ^j. The phonon fre- 
quencies at the r and A points agree remarkably well. 
Some differences appear at symmetry points K and M. 
Such diffemnces are also present between the results of 
Alfe et al£3 and those obtained by Soderland et aZr— us- 
ing a generalized pseudopotential parameterizationE3 of 
FP-LMTO calculations. The differences between our LR 
results and those of Alfe et al£3 are not large enough 
to cause significant differences in thermal properties and 
electron-phonon coupling. The smooth solid lines in Fig. 
H correspond to spline fits to the calculated frequencies 
(solid circles). Due to the small number of calculated fre- 
quencies the shapes of the lines presumably representing 
the bands at the zone boundaries could be incorrect. The 
connections of the calculated points with lines and band 
crossings in Fig. ^ were determined by examining the 
phonon eigenvectors. However, the number of q-points 
considered along each symmetry direction was at most 
four and often three or less. No intermediate q-points 
along the K - M and L - H were among the mesh of q- 
points for which the dynamical matrix was calculated. 
Thus the possibility of errors in band crossing cannot be 
ruled out. 

The dispersion curves at various pressures are similar, 
except for an overall scale factor, essentially representing 
the gradual broadening of the bands with increasing pres- 
sure. This is reflected in Fig. || and also in the phonon 
density of states for various lattice parameters shown in 
part (b) of Fig. |^. For the smallest lattice parameter 
considered by us the upper band edge lies around 670 
cm-i or 20 THz (Fig. | and Fig. | (b)). 

We have computed both the Eliashberg spectral func- 
tion 
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FIG. 2: Phonon frequencies of hep Fe calculated via FP- 
LMTO-LR method for two different lattice parameters: (a) 
4.0 and (b) 4.4 a.u., and the ideal c/a ratio, ^8/3. The solid 
circles denote the calculated frequencies and the solid lines 
represent spline fits through these calculated values. 

and the transport Eliashberg functiorHS 
1 



(7) 



where the subscript F denotes the Fermi surface, the 
angular brackets denote the Fermi surface average, vp 
denotes the Fermi surface velocity, g^'^, is the electron- 
phonon matrix element, with v being the phonon polar- 
ization index and k, k' representing electron wave vectors 
with band indices i, and j, respectively. 

For most of the lattice parameters considered by us 
the Eliashberg spectral function a^F and the transport 
Eliashberg function a^^F both follow the same frequency 
variation as the phonon density of states. In Fig|| we 
show the phonon density of states and the two Eliashberg 
functions together with the phonon dispersions for the 
lattice parameter 4.6 a.u., close to the equilibrium value 
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of 4.615 a.u. 

Some deviations in the frequency dependence of the 
a^F function from that of the phonon density of states 
appear at higher pressure. The deviation is most pro- 
nounced between the lattice parameters 4.4 and 4.2 a.u. 
in our calculation. In Fig.^we show the Eliashberg spec- 
tral functions and the phonon density of states for three 
different lattice parameters. The peaks in the calcu- 
lated phonon density of states at ambient pressure (lat- 
tice parameter ~ 4.6 a.u. in our calculations) are at 
190 cm-i (24 meV) and 315 cm"! (39 meV), and these 
agree very well with the results from neutron scattering 
experimental as well as with recently reported results, 
obtained from the pleasured energy spectra of inelastic 
nuclear absorption.Ej The peak positions in the calcu- 
lated results for higher pressures are at somewhat lower 
frequencies (by about 5 meV, which is within the exper- 
imental resolution) than-those from the inelastic nuclear 
absorption experiment.^ However, such differences be- 
tween the calculated and measured frequencies are com- 
mon, given the difference between the experimental and 
theoretical values of the lattice parameters at various 
pressures. 

The Hopfield parameter rj {N{0){l'^)), which is the 
electronic part of the electron-phonon coupling, shows 
an above average increase between the lattice param- 
eters 4.5 and 4.4 a.u., due to an increased coupling 
for the longitudinal acoustic phonons with wave vectors 
around the middle of the F- K symmetry line. This in- 
creased (above average) coupling is found to persist up 
to at least 4.2 a.u., but diminishes to normal (average) 
value around the lattice parameter of 4.0 a.u., where the 
pressure-stiffening of the lattice vibrations reduces the 
electron-phonon coupling and superconductivity disap- 
pears. The general trend is as follows: the Hopfield pa- 
rameter ry grows steadily with increasing pressure with a 
rapid change between the lattice parameters 4.5 and 4.4 
a.u. The phonon frequencies move upward with increas- 
ing pressure, with no phonon branches showing any soft- 
ening. However, between lattice parameters 4.6 a.u. and 
4.4 a.u. (perhaps 4.3 a.u.) the increase in the Hopfield 
parameter dominates the change in the electron-phonon 
coupling parameter Xph = r]/M{uj^). In this range Xph 
increases despite a decrease in A'^(O) and an increase in 
(w^). Below 4.3 a.u. lattice vibrations stiffen rapidly, 
lowering the value of Xph- In Table |l| we summarize our 
results for the pressure-dependence of the phonon prop- 
erties and electron-phonon coupling. Since the Hopfield 
parameter 77 is often calculated using the rigidjauffin-tin 
(RMT) approximation of Gaspari and Gyorffy ,e2I in Table 
y we have also presented the results for rj obtained via 
the LMTO-ASA implementation of RMT^rigid atomic 
sphere or RSA) as given by Glotzel et alE3 and Skriver 
and Mertig.eHl These values are in agreement with those 
given by Mazin et-a/.Jij but differ significantly from the 
results of Jarlborgn3 (judging from the quoted values of 
Xph and the Debye frequencies). Our results indicate that 
depending on the lattice parameter the RMT/RSA ap- 



proximation underestimates the Hopfield parameter by 
15-45%. Also the variation of rj with lattice parameter in 
the RMT approximation is much smoother than in the 
linear response calculation, as the former fails to capture 
the above average increase around the lattice parameter 
4.4 a.u. In Table |ll| we have presented the results for lat- 
tice parameter 4.7 a.u. merely for comparison with other 
lattice parameters, and not for comparison with exper- 
iment. The strong electron-phonon coupling (stronger 
than that at a =4.6 a.u.) is of no experimental conse- 
quence since, (i) at this lattice parameter the system is at 
a negative pressure, not accessed by experiment; and (ii) 
our theoretical calculations show that at this expanded 
volume the system is most likely antiferromagnetic. 



IV. CRITICAL TEMPERATURE 



A. General relations 

The linearized Eliashberg equations at the supercon- 
ducting transition temperature Tc of an isotropic system 
are (see, e.g., Ref.Ell): 



Z(zc^„) = 1 + — Vl^+(n-n')sign(n'),(8) 

n 



wherein = 7^^c(2f^^- 1) is a Matsubara frequency, A(ia;„) 
is an order parameter and Z{iujn) is a renormalization 
factor. Interactions 

W+[n - n') = Xph{n - n') + Xsf{n - n') + Snn'iln + 7m), 



and 

W-{n - n') = Xph{n - n') - Xsf{n - n') - ^j.*{uJc) + Snn'iln - 7m), 
contain a phonon contribution 



Xph{n - n ) 







(cj„ - UJn'Y + 



where a^{uj)F{Ld) is the so-called Eliashberg spectral 
function, and a contribution connected with spin fluc- 
tuations 



Xsf{n - n) 



00 ^,,2 



djjj^P{u) 



P(w) is the spectral function of spin fluctuations, related 
to the imaginary part of the transversal spin susceptibil- 
ity X±{^) as 



P{oj) = (|5kk'|^Imx±(k,k',Lj) 
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States/cm 



FIG. 3: Phonon spectrum, density of states and the Elias hber g spectral function F and the transport Ehashberg function 
a^'^^ F for hep Fe at the lattice parameter 4.6 a.u. (c/a=^8/3). The equilibrium (minimum energy) lattice parameter is 4.615 
a.u. 



and ii*(luc) is the screened Coulomb interaction 



1 + fi\n{E/uJc)' 



(9) 



where /i = {N{0)Vc) pg is the Fermi surface averaged 
Coulomb interaction. E is a characteristic electron en- 
ergy, LUc is a cut-ofF frequency, usually chosen ten times 
the maximum phonon frequency: lOc — IOlu^j^^. 7„ = 
l/2T„,7m = 1/2t,„ are scattering rates for nonmagnetic 
and magnetic impurities, respectively. 



B. Phonons only 

In order to compute Tc we use the calculated Eliash- 
berg spectral function along with the following proce- 
dure to determine the Coulomb pseudopotential ^*{u!c)- 
We start by assuming fi — 0.5. A value greater than 
0.5 for /i would lead to magnetic instability (see, e.g. 
Ref.EsI). With E = ujpi, the electron plasma frequency 



(see RefJH), 



0.5 



1 + 0.5\n{ujpi/u!c) 



Thus from the calculated phonon frequencies and plasma 
frequencies we obtain /x* for all lattice parameters, with 
the cut-off frequency lOc assumed to be ten times the 
maximum phonon frequency. This procedure gives us 
the maximum possible values of fi*(LUc). 

One of the most widely used, expressions for Tc is 
given by the Allen-Dynes formo of the McMillan for- 
mula (Eq.|l|), where the prefactor 8d/1.45 is replaced by 

ph 



-71.2. 



\ph = ^ / dcjci {uj)F{uj)/lu 

10 



is the electron-phonon coupling constant, wJJ^" is a loga- 
rithmically averaged characteristic phonon frequency 

uj^^^ — exp < - — f dwa'^ {lj)F{lj) In uj / Lj \ , 
I \h Jo J 
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TABLE 11: Hopfield parameters from the linear response calculation rj and the rigid muffin-tin (atomic sphere) approximation 
f] (RMT/RAS), mean square electron-ion matrix element {I^), calculated average plasma frequencies ujpi, logarithmic average 
phonon frequencies win, cutoff frequencies ujc, Coulomb pseudopentials for Eliashberg equation {jj,*{ujc)) and McMillan formula 
iHin); electron-phonon coupling parameters Aep, calculated critical temperatures {T^^^'^) and superconducting gaps (Ao) from 
the solution of the Eliashberg equations and the critical temperatures from the McMillan formula (^)(rj^'^^) for various 
lattice parameters a. 
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FIG. 4: Phonon density of states and the Eliashberg function 
for hep Fe for three different lattice parameters (c/a=^8/3). 



and 



1 + ^*(wc)ln(a;c/c^fn ) 



is the Coulomb pseudopotential at this frequency. Our 
calculations show that for different plasma frequencies 
and characteristic phonon frequencies ji* for all lattice 
parameters lies in the range 0.13 — 0.14, which is typical 
of conventional superconductors. 

In Fig.^ we show the transition temperatures calcu- 
lated as a function of volume per atom using Eliashberg 



equations and the McMillan formula. The effects of fer- 
romagnetic and antiferromagnetic spin fluctuations, dis- 
cussed in the following subsection, are also shown via 
three additional curves. The symbols denote the calcu- 
lated values of Tc and the lines are spline fits through the 
calculated values. 



C. Contribution from spin fluctuations 



Superconducting transition temperatures calculated in 
the previous subsection are based on the maximum possi- 
ble estimates of the Coulomb pseudopotential /i*. Thus 
Tc, based on s-wave electron-phonon interaction only, 
cannot be less than 4.5 K, and a value as high as 7- 
8 K is reasonable according to the linear response re- 
sults. The highest transition temperature obtained in 
the experiments is 2 K. A more important difference be- 
tween the calculated and the experimental results is the 
range of volume/pressure over which superconductivity 
appears. The calculated range is much broader than the 
experimental one (see Fig.^). It is then natural to explore 
the effects of spin fluctuations on both, the magnitude of 
Tc and the pressure/volume range of the superconduct- 
ing phase. Since the calculation of the spin succeptibil- 
ity is rather complicated, we restrict ourselves to simple 
models of ferromagnetic and antiferromagnetic spin fluc- 
tuaticps for an isotropic system, as proposed by Mazin 
et a/..Ej In a T— matrix approximation for the uniform 
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55 65 
Volume per atom (bohr ) 

FIG. 5: Calculated transition temperatures for various vol- 
umes per atom. The experimental results are also shown. The 
experimental pressure versus Tc results were transferred injfeei 
volume versus Tr. results using the data from Mazin et aZtfl 
(see also Refs.QH). The legends SPFl, SPF2, and AFMSF 
are described in section IV. C. Dashed vertical lines show re- 
gions, where hep Fe is known to be antiferromagnetic (AFM) , 
and where it is either nonmagnetic (NM) or antiferromagnetic 
(AFM). 



electron gas one can obtain the relation (see Refs.p? 48) 



P{uo) = iV(0) / rfg^ j -ilmx±(<Z,c^) 



where 



-ilmx±(g,t^) 

(l - /iV(0) - IN{Q) 



An integration of P(a;) (see section A.) leads to the spin 
fluctuation coupling parameter 



A 



= aiV(0)/hi — 



1 



N{Q)I' 

where the constant a is of order unity. One can define 



(10) 



exp - — / duiP{u})\rLU} /lo 
^sf Jq 



(11) 



as a characteristic spin fluctuation frequency, which 
should vanish near the magnetic phase transition, pp, 
and vp are the Fermi momentum and velocity, respec- 
tively. The product pfVf can be replaced by 2Ep and 
estimated from the location of the Fermi energy with re- 
spect to the bottom of the band. 
If we resort to the approximation 

P(w) = {XsfUJsf/2)S{uj - UJsf) , 



then for ajgf ^ u!ph we obtain an extension of the McMil- 
lan formula, similar to the one used in Rcf.[ll 



ph 

1.2 



exp 



1.04(1 + Xph + Xsf) 



Xph - Xsf - /^*[1 + 0.62(Ap/i + Xsf)] 



(12) 



In reality the spectrum P{lu) is distributed from zero 
up to electronic energies. Near the phase transition the 
characteristic frequency is comparable to the characteris- 
tic phonon frequencies. An appropriate treatment of the 
broadness of the spectrum P{uj) leads to the replacement 
of the ojf^'^ in the above expression by 



with the exponent v (see Ref.H8) given by 



X2 



(13) 



(A,, - Xsf) [XpH - Xsf + ,X';if i^KnV'^./] 



ph 



In the uniform electron gas approximation the constant 
a in Eq. (|l^) is of the order of unity. But such a high 
value of a in our calculation would cause the critical tem- 
peratures to vanish for all lattice parameters {Xph < Xsf). 
It is evident that the uniform electron gas approximation 
would be inappropriate for a transition metal like iron. 
Hence we use the following approach: we consider a as 
a fitting parameter to get a Tc =2 K, the experimen- 
tal value, for the lattice parameter 4.4 a.u. (volume per 
atom ~ 60 a.u. The two sets of Tc versus lattice parame- 
ter results obtained this way are shown in Fig. (^ and are 
labeled as SPFl and SPF2, respectively (the lowermost 
curves). In particular, SPFl refers to the case where 
Eq. ^ is used with Xsf given by Eq. dl^); and SPF2 
refers to the case where oj from Eq. ( |13|) r eplaces uif^ 
in Eq.(|l|), with A,/ still given by Eq. (|o|). In Table 
III, the spin fluctuation coupling parameters associated 
with the results SPFl and SPF2 in Fig. | are labeled 
as As/1 and As/2, respectively. The values of the param- 
eter a (see Eq.(|l^)) for the two cases SPFl and SPF2 
are 0.101 and 0.029, respectively. Calculated transition 
temperatures for the two models SPFl and SPF2 are de- 
noted by Tj/ and T*/ in Table III , which also shows the 



characteristic spin fluctuation frequencies us^^^ for various 
lattice parameters. 

For antiferromagnetic spin fluctuations the spin sus- 
ceptibility has a maximum at q — > Q*, and the averaging 
over the Fermi surface leads to 



A 



A(q^Q*)/ 



i-x(q^Q*)/' 

If, according to Mazin et we suppose x(q Q*) 
bN{0), then 
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TABLE III: Spin fluctuation effects: ferromagnetic spin fluctuation coupling parameters \sfi,Xsf2 (see text in section IV. C. 
for details); antiferromagnetic spin fluctuation coupling parameter X'^j, the characteristic spin fluctuation frequency ti^f/, and 
the corresponding critical temperatures ,T^2 ^ ^^nd T^{f (see text in section IV. C. for details). 
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j_ a'bN{Q)I 

\-hN{o)r 



(14) 



Parameter b can be estimated from the condition of the 
antiferromagnetic instabihty bN{Q)I 1. This leads to 
6 ^ 1.5 which is close to the value in RefjlJ. Taking 
this value and using a' as a fitting parameter we obtain 
the result plotted in Fig. ||. a' = 0.0032 reduces the 
maximum Tc (at 4.4 a.u.) to the experimental value, 
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), with Lu^^ replaced 

af 



2 K (for simplicity we used Eq.( 

by u! given by Eqns. ( pT| ) and (13]), and A"^ replacing 
As/). The corresponding results are plotted in Fig.(||) 
and are labeled as AFMSF. The values of Tc and spin 
fluctuation coupling parameters are also shown in Table 
III, labeled as T^^ and A°J, respectively. It is evident 
that the volume dependence of Tc obtained this way is 
very similar for ferromagnetic and antiferromagnetic spin 
fluctuations. 



c o Eliashberg equations 

n a AG theory 




FIG. 6: Variation of the transition temperature with the 
scattering rate of magnetic impurities l/2rm for lattice pa- 
rameter a =4.4 a.u. The legend AG stands for solution of the 
Abrikosov-Gor'kov expression given by Eq.(|l5|). 



D. Magnetic impurities 

A lowering of the critical temperature could also be 
caused by the presence of magnetic impurities. It is well- 
known that the nonmagnetic impurities canceLout from 
the Eliashberg equations (Anderson theorerrCJ )p-,whilc 
the magnetic ones lead the pair-breaking effectsrS The 
central idea is that near a magnetic transition spin- or- 
dered clusters appear, and these can scatter electrons 
very effectively. We have calculated the effect of such 
impurities on the critical temperature for the lattice pa- 
rameter a = 4.4 a.u. (see, Fig. |^) by considering various 
different scattering rates l/2rm in the Ehashberg Eq.(||). 
For comparison we have also calculated the change in Tc 
by using the renormalized Abrikosov-Gor'kov (AG) ex- 
pression (see e.g., section 15 in Rcf.^ 

ln(r,o/T,) ~ iP [1/2 + (l/2T,„)^Te(l + A^p)] - ^(1/2), 

(15) 



where Tco is the critical temperature without magnetic 
impurities, i/jix) is the digamma function and 'ip{l/2) is 
related to the Euler constant 7 as -0(1/2) — —7 — 2 In 2. 
The difference between the Eliashberg and the AG re- 
sults is due to the rather broad phonon spectrum which 
necessitates appropriate treatment of strong coupling ef- 
fects. 

In order to reduce the critical temperature at the 
lattice parameter a = 4.4 a.u. to the experimental 
value of 2 K it is sufhcient to assume a scattering rate 
I/tto w 3.0 cm^^. With the calculated average Fermi 
velocity vp — 2.547 x 10^ cm/sec, this yields a mean free 
path I ~ 28.3 X 10~^ cm. The closer the magnetic in- 
stability, the larger is the probability (rate) of magnetic 
scattering, which leads to more enhanced suppression of 
Tc. 
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E. 



p-wave pairing 



Magnetic ordering (as well as an external magnetic 
field) favors the triplet p-wave pairing, similar to that 
found in superfluid '^He. In order to estimate Tc for p- 
wave pairing we adopt the following simplified approach. 
We consider the extension of Eq. (||) for the Z*'* spherical 
harmonic channelEll: 



1.2 



exp ■ 



Y^^cxp. 



1 + Ap/i + As/ 
g{Xph + Kf) . 



(17) 



Z(iu^) = 1 + — y Ty|°^(n-n')sign(n'), (16) 

^ Z{lUJn')UJn>\ 
n' 

where d*-'-* for / = 1 is the p-wave order parameter, and 



where uj is given by Eq. {[^, and we have used the 
relation: A^j = gX^j: = gXsf- Note that this equat 
tion is the same as equation (3.6) of Fay and Appel,Ej 
except that the term Ap^ is^ absent from the exponent 
in their expression for A small value of the pa- 

rameter g and a rather strong phonon contribution to 
the numerator in the exponent in Eq.(p^ lead to small 
values of Tc ^ 10~^K for the p— wave pairing in con- 
trast to the conclusion reached in Ref.^. The value 
Tc ^ IQ— similar to that obtained by Allen and 
MitrovicEi for p-wave superconductivity in Pd. If the 
m/^'Vti — \' ^ \' a i"^, _u^, \assumption M^(^) = gW^°^ is valid, then the inclusion 

01 antiierromagnetic spin nuctuations (replacing As/ by 
As/) would lead to similar results. 



W^\n - n') = A^,,(7i - n') + {-lYxifin - n') + 5ioX , 



SUMMARY AND CONCLUSIONS 



where X = — H* {u>c) + Snn' {in —7m)- The kernel W with 
a general index I is defined as the Fermi surface (FS) 
average 



,n — n )dW 



of the harmonic of the momentum-dependent interac- 
tion WCk, k', n — n), while 



W'-°^ = (M^(k,k', 



k,k'e_F5 ' 



We assume that the Coulomb interaction and impu- 
rity scattering are isotropic. The simplest approxima- 
tion then is to use ^ gW'-°\ where the param- 
eter g describes the anisotropy of the interaction (see, 
e.g., Refs.^ 52). A difference of the factor g from unity 
leads to strong pair-breaking effects. In general (see., e.g. 
Ref.lll), for Z = 1 d^ ' ~ Vi?, the first odd Fermi surface 

harmonicS In this case g = ^ph/Aph = ^ph^ /\h (see 
notations in Refs.^ 20). The phonon constant Ap™-* is 
relevant to the transport Boltzmann equation (see, e.g., 
Ref.EOh. From the linear response calculation we obtain 
g — Xpf^^'/Xph for all lattice parameters. For a=4.4 a.u., 
we obtain g=0.238. 

With the assumption W'^^^ — gW'^'^^ it is possible to 
estimate Tc from a McMillan-like formula. An expres- 
sion for the critical temperature can-be written in a form 
similar to that used by Mazin et aZJiJ: 



The state of hep Fe under moderately high pressure 
(< 60 GPa) is currently riddled with controversial re- 
sults that need to be resolved and understood. The ini- 
tial experimental result by Shimizu et alu indicating su- 
perconcLuctivity between 15 and 30 GPa has been verified 
recentljfl for a pressure of 22.5 GPa, where the maximum 
Tc of 2 K was observed. Experimental measurements of 
Raman spectra are also reported. Ej'LZI The observation 
of a second Raman peak besides that due the Raman 
active i?2g mode in hep Fe has bee n assigned to disor- 
der induced phonon scattering ,E3EZl although, based on 
density functional calculation, the possibility of antifer- 
romagnetic order up to a pressure of approximately 60 
GPa has also been suggested.Ea Normal state resistivity 
measurements by Jaccard et al. show a T^^^ tempera- 
ture dependence at low temperature, which is predicted 
by a nearly ferromagnetic Fermi-liquid model. The ear- 
lier Mossbauer study of hep Fe under pressure had failed 
to reveal any local magnetic moment ,c3'Ej and results of 
theoretical calculations are dependent on the treatment 
of the exchange and correlation potentials used in the 
calculations. In view of thia. we have adopted the same 
approach as Mazin et a/Jlil: we assume nonmagnetic 
state under pressure and examine the electron-phonon 
coupling as a function of the lattice parameter. 

The results of our study can be summarized as follows: 
(i) The Hopfield parameter i] increases steadily with pres- 
sure, showing a wider variation for the linear response 
calculation than obtained via the RMT approximation 
(both in our J,MTO-RMT calculation and in that of 
Mazin et alS^ ). (ii) Below volumes '^50 a.u. per 
atom (above estimated pressures ~ 160 GPa) phonons 
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stiffen rapidly, bringing the Tc down (somewhat faster 
than what is suggested in Ref.^^. (iii) T^'s based on the 
s-wave electron-phonon couphng, and maximum possi- 
ble estimates of /i* are higher than the experimental val- 
ues (iv) The range of volume where superconductivity 
appears is much broader in the calculations than what 
is observed, in agreement with the result of Ref, 14. (v) 



Inclusion of ferromagnetic/antiferromagnetic spin fluctu- 
ations, and scattering from magnetic impurities can all 
bring the calculated values of down to the range of 
observed values, but cannot substantially improve the 
agreement between the calculated and the experimen- 
tal pressure/volume range of the superconducting phase 
(Fig.^). (vi) A simplified treatment of p-wave pairing 
due to electron-phonon and spin fluctuation interactions 
yields a very small Tc (< 0.01 K), in contrast with the 
claim made in Ref.|2^. 

The role of impurities remains somewhat puzzling as 
well as of vital iuterest at present. The initial results 
of Shimizu et alh showed very high normal state resid- 
ual resistivity AQfiVLcva), indicating the presence of 
substantial de&cts in the samples studied. More recent 
measurementg3 with purer samples (with residual resis- 
tivity 50 times smaller) show the same maximum at 
the same pressure. If the impurities in the earlier stud- 
ied samples were magnetic this would rule out electron- 
phonon s-wave coupling as the primary mechanism of 



superconductivity. However, Jaccard et a/.BI also notice 
that superconductivity in their samples is unusually sen- 
sitive to disorder, developing only when the electronic 
mean free path exceeds a threshold value. They find the 
normal state resistivity to be characteristic of a nearly 
ferromagnetic metal, but our calculation of p-wave su- 
perconducting Tc in the presence of ferromagnetic spin 
fluctuations and electron-phonon interactions yields val- 
ues less than 1Q~^ K, similar to that for Pd obtained by 
Allen and Mitrovico 

Our calculations seem to rule out both s- and p-wave 
superconductivity in hep Fe. Electron-phonon mediated 
s-wave superconductivity should persist, in severe dis- 
agreement with experiment, beyond pressures of 200 GPa 
even in the presence of spin fluctuations. The possibil- 
ity of d-wave superconductivity mediated by antiferro- 
magnetic spin fluctuations remains to be explored. A 
d-wave superconductivity would be consistent with the 
observationcl that superconductivity in hep Fe seems to 
be extreemly sensitive to disorder. Low temperature spe- 
cific heat measurements can further clarify this issue. 
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